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Abstract. The properties of spin polarized pure neutron matter and symmetric 
nuclear matter are studied using the finite range simple effective interaction, upon 
its parametrization revisited. Out of the total twelve parameters involved, we now 
determine ten of them from nuclear matter, against the nine parameters in our earlier 
calculation, as required in order to have predictions in both spin polarized nuclear 
matter and finite nuclei in unique manner being free from uncertainty found using the 
earlier parametrization. The information on the effective mass splitting in polarized 
neutron matter of the microscopic calculations is used to constrain the one more 
parameter, that was earlier determined from finite nucleus, and in doing so the quality 
of the description of finite nuclei is not compromised. The interaction with the new set 
of parameters is used to study the possibilities of ferromagnetic and antiferromagnetic 
transitions in completely polarized symmetric nuclear matter. Emphasis is given to 
analyze the results analytically, as far as possible, to elucidate the role of the interaction 
parameters involved in the predictions. 
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1. Introduction 

The study of nuclear matter and hnite nuclei in a given model is a subject of 
contemporary interest in the area of nuclear research. The most fundamental ab initio 
calculations of Dirac-Brueckner-Hartree-Fock (DBHF), Brueckner-Hartree-Fock (BHF) 
and variational types [Tl|2l|3llll[5l|6l|7llHl[9l UHl [HI [I2] start from a Hamiltonian 
which is adjusted to reproduce the nucleon-nucleon (N-N) scattering phase shifts and 
properties of few nucleon bound systems. The predictions in the regime of nuclear 
matter (NM) of this kind of ab initio calculations are usually considered as a standard. 
However, the extension to hnite nuclei of microscopic calculations has severe constraints 
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due to the much involved theoretical and computational procedures. Mean held model 
calculations [HI H [ISmsi [H] are very popular to deal with hnite nuclei for their 
relatively simpler computational requirements and analytical advantages as compared 
to ab initio calculations. 

Within such a kind of models the Relativistic Mean Field (RMF) model, which 
uses an effective Lagrangian constructed by considering meson exchange, enjoys an 
advantageous position for its successful application to structure and reaction studies of 
hnite nuclei including stable as well as superheavy nuclei [ISl UHl UHl EQl EU [22], l2^ l2l]. 
In the non-relativistic domain, the Skyrme interactions [TH |25l |26l [23 [28] [29] deserve a 
similar status for their wide applications to hnite nuclei calculations due to the analytical 
simplicity because of the zero-range of these forces. A considerable progress has also 
been made in the non-relativistic domain to develop mean held models using hnite 
range interactions. The Gogny (SO] El] EH [33], M3Y [SU [35] and the hnite range simple 
ehective interaction (SEI) [TG] [36] E3 forces are examples of this type of interactions. 
In addition, it is worth mentioning the recent development of hnite range functionals 
using the ehective held theory methodology applied to low-energy nuclear physics [38] . 

The equation of state (EOS) and the momentum dependence of the mean held are 
two important aspects in the studies of NM (a compilation of several acronyms used in 
this paper is provided in Table 1). However, ehective mean held models may predict 
results in NM that do not necessarily agree with the results provided by microscopic 
calculations. For example, microscopic calculations predict a larger ehective mass for 
neutrons than for protons in neutron-rich isospin-asymmetric nuclear matter (ANM) 
[61 EU SQl E], which is also the currently accepted point of view, as substantiated 
by the experimental results of the energy dependence of the Lane potential [441 [45] . 
However, mean held calculations based on ehective interactions in both the relativistic 
and the non-relativistic domains do not always fulhll this trend [lUialSniiT]. This 
points out to the fact that the momentum dependence of the mean held in NM has not 
been taken as a prerequisite in the htting of the parameters of the RMF and most of 
the ehective non-relativistic models. As a consequence, there is no reason a priori why 
these models should follow the trend of momentum dependent properties, such as mean 
held, ehective mass splitting, etc. as provided by microscopic calculations in NM. The 
parameters of ehective relativistic and non-relativistic models are usually constrained 
from the empirical values of some NM properties, such as energy per particle e(p), 
incompressibility K{p), symmetry energy Es{p), etc., computed at saturation density 
p=Po and to some selected experimental data of binding energies (BE) and radii over the 
periodic table. The empirical values of po, e(po) and Eg^po) used in these models vary 
within the ranges 0.17±0.03 fm“^, —16±0.2 MeV and 33±5 MeV, respectively. With 
htting protocols of this type, the momentum dependence of the mean held is completely 
predicted by the model and can or cannot reproduce the tendency exhibited by the 
microscopic calculations. 

The momentum dependence of the mean held is a fundamental property [49l [50] 
and it should not be left open to assume an arbitrary behaviour. This momentum 
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Table 1. Several acronyms used in the text. 


Acronym 

Meaning 

SEI 

Simple ehective interaction 

NM 

Nuclear matter 

BE 

Binding energy 

SNM 

Symmetric nuclear matter 

ANM 

Asymmetric nuclear matter 

PNM 

Pure neutron matter 

PPNM 

Polarized pure neutron matter 

CPNM 

Completely polarized pure neutron matter 

EM 

Ferromagnetic 

AFM 

Ant iferro magnet ic 

CSNM 

Ferromagnetic completely polarized SNM 

CASNM 

Antiferromagnetic completely polarized SNM 


dependence of the mean field, as extracted from the analysis of nncleon-nncleus 
scattering data [501 EB ESI EH] is explicitly taken into account in the htting procedure 
of the parameters of the hnite range simple effective interaction (SEI) which has been 
used in NM studies [17l[56l[57l EH EH] . In a recent work [37], the SEI has been extended 
for studying ground-state hnite nuclei properties. It should be pointed out that the 
SEI depends on eleven parameters (apart from a spin-orbit strength parameter in the 
case of hnite nuclei) of which nine could be obtained from the studies of ANM. The 
parameters responsible for the momentum dependence of the nucleonic mean helds in 
ANM are decided by utilizing appropriate constraints with care that these predictions 
are not changed while deciding the rest of the parameters. The two open parameters of 
SEI, to and Xq, and the spin-orbit strength Wq are left to reproduce a few magic nuclei. 
Within this protocol to determine the parameters of the SEI, the microscopic trends of 
the properties in ANM could be reproduced and, at the same time, the known binding 
energies and charge radii of even-even spherical nuclei could be described with a quality 
similar to other successful mean field models of relativistic or non-relativistic type (see 
table 3 of [6^) [37] • However, while examining the results in spin polarized NM it is 
found that the determination of the parameters to, xq and Wo from hnite nuclei does 
not yield unique predictions in spin polarized NM. In this work our objective is to study 
the spin polarized NM using the SEI in consonance with the microscopic predictions to 
remove the aforementioned ambiguity. This is done by using the microscopic information 
on the momentum dependence of the mean held in spin polarized pure neutron matter 
(PPNM) to hx the xq parameter, leaving R and Wo only to be determined exclusively 
from hnite nuclei. Our procedure enables to constrain the strengths of the interaction in 
the four basic channels of the N — N interaction, namely, the singlet-even triplet- 
even triplet-odd and singlet-odd channels as well as the predictions 
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in spin polarized NM, in a unique manner for a given EOS. In this new method of 
determining the parameters, the earlier predictions in ANM do not change and the 
hnite nuclei results are also reproduced within reasonable accuracy. Thus the obtained 
SEI can be used in the study of both isospin and spin polarized NMs as well as in 
hnite nuclei. In section 2 we outline the formulation of spin polarisation in symmetric 
nuclear matter (SNM) and pure neutron matter (PNM) using the SEI. In this section 
we also discuss the determination of the parameters from ANM and hnite nuclei, as well 
as the underlying uncertainty that manifests in the predictions in the spin channel. In 
section 3 the procedure adopted to remove the uncertainty is worked out and the new 
SEI parameter sets for EOSs having different NM incompressibility are obtained. The 
predictions in the spin channel in PNM and SNM are discussed and compared with the 
results of other mean held models as well as with microscopic calculations. Section 4 
contains a brief summary and conclusions of our analysis. 


2. Formalism 


The hnite range simple effective interaction (SEI) used in the present work is given as 


'^effir) = to{l + XoP^)5{r) 

(IF + BP^ - HPr - MP^Pr) f{r), 


+ ?(1+X3P.) 


5{r) 


( 1 ) 


where, /(r) is the functional form of the hnite range interaction containing the single 
range parameter a and is taken to be of Gaussian form, The other terms have 

their usual meaning. The SEI in equation ([1]) has 11 parameters, namely, b, to, xq, 
G, X3, 7 , a, W, B, H and M (the spin-orbit strength parameter Wq will enter in the 
formulation of hnite nuclei). The complete study of ANM requires the knowledge of 
altogether nine parameters, namely, b, 7 , a, Eq and Eq^ with, the indices 

” 1 ” and ”ul” depicting interactions between pairs of isospin-like and unlike nucleons, 
respectively [571 [58] . For the sake of simplicity, the formulation has been based on the 
fact that the range between a pair of isospin-like or unlike nucleons is the same but 
they differ in their strengths. The connection between the parameters of ANM and the 
interaction parameters is given in the earlier works [SUET]. Here we shall write them 
in terms of the strengths of the hnite range part of the N-N interaction in the four 
states, namely, singlet-even Vq^, triplet-even Vq^, triplet-odd Vq^ and singlet-odd 
instead of W, B, H and M for the sake of convenience of the discussions in this 
work. They read as, 

el' = (2 + 1 „) + ^ {Veie + 31^® + 3V™ + 1®®°) j f(,r)d\ (2) 

= jPo (1 - xo) + ^ + 317°) j f(.r)d‘r (3) 


( 4 ) 
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£ 


I 

7 


- xs) 


(5) 


£ 


ul 

ex 


f (7/'= + 3Vr - 3Vr - Vf°) 


f(r)d^r 


( 6 ) 


4 = J (V,f‘' - 3Vi™) / (7) 

The knowledge of the nine parameters that characterize the ANM can be obtained from 
the independent studies of spin saturated SNM and PNM. 

In SNM the numbers of neutrons and protons are equal. In the case of completely 
polarized SNM, there are two possibilities for the spin polarisation. One of them 
corresponds to the situation where the spins of neutrons and protons are aligned in the 
same direction (CSNM), referred to as ferromagnetic (FM) type. The other possibility 
corresponds to the situation where the spins of neutrons and protons are oriented in 
opposite directions (CASNM), referred to as anti-ferromagnetic (AFM) type. For the 
SEI, the energy per particle (obtained as the ratio of the energy density H with the 
density p) in SNM, CSNM, and CASNM is given, respectively, by 


( 8 ) 


(9) 


e(p) = 


H{p) 


+ 


~ ~ lOM 

\^ex ' ^exJ 


+ 


(4 + ^0 ) 


4po 


P + 


(4 + O 


4po 


pj{kf) 


P 


P(R) Y 

1 + bpiK)) 


-'pol 


{p) = 


hUp) 

p 

-/s 






^Is 


lOM 


+ 7T^P + 


2 po 2 p| 


7+1 


P 


+ ) 


P(R) Y 

i + MR)y 


where 



and A=4 The expressions of J{ki) in equations (jH]), (jO]) and flTOll can be obtained from 
equation (ITT]) by using ki = kf = (|vr^p)3 and ki = = (37r^p)3, where kf is the 

Fermi momentum in SNM and is the Fermi momentum in CSNM and CASNM. 
The new parameters appearing in equations ([9]) and fiTOj) with index “Is” and “las” are 
given as. 


£ 


Is _ 

0 “ 


( 1 +^o) + ^ 



( 12 ) 
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4 ' = + ^3) (13) 

4 '. = f (vr - 314 °) /( 14 ) 
4" = |l>» + f (4"'" + + sv;™ + Vi*®) j f{r}£r ( 15 ) 

4" = (16) 


_/a5 


“ ^ V'^o 


+ - 515 /° _ v/°) 


f(r)d^ 


(17) 


Similarly, the energy per particle in PNM, PPNM and completely polarized PNM 
(CPNM) can be given for the SEI as, 




(p) = 


+ 


H^jp) 


3ffkl Eq 

- - H- -p 

lOM 2po 

l + 6p(R)) 


p(R) 


( 18 ) 


and 


^ o \ _ ^pol(Pn‘>P^ P'n-d) _ 1 r ^^^{KiuPnu + k‘^dPnd) ,2 , 5 

^polyPiP^p) in/i/f "*^0 yrnu' rnd) 

^ p pi lUM zpo 


+ 


^l,ul 

-0 


, f P r 2 , 2 5 ^ 7 “^ ^ P(R') ^ 

PnnPnd + [^^{Pnu + Pnd) + -^PnuPnd) ^ J 


PO 2 pg Pq 

+ + PndJ^^nd)^ 

^nd pkjiu 

2 , 


+ 




l,ul 


4povr2 


Pn 


I{k,knu)k dk + Pnd / I{k,knd)k dk 


'0 


'0 


( 19 ) 


^cpnm (P) = 


KnmiP) _ , £0 


P 


lOM ^ 2po^ 


-Z,Z 


p(R) \ 


2 p 4 ’ V1 + 6 p(R) y 2 po 


’ + dfpJ(P/'), 


( 20 ) 


where, in equation (IT^ . 


I{k, hi) = 


3 A 3 rA 




A: ^ 


' k+kj ( k — kj \2 

) — e~v“^) 


+ v/^ erf 


k + h 
A 


— erf 


k — ki 
A 


( 21 ) 


for /cj = knu, knd- One has /c„=(37r^p)3 for the Fermi momentum in PNM, and 
knu(nd) = {Q'^‘^Pnu{nd))^ and A:P°*= (67r^p)5 for the Fermi momentum in PPNM and CPNM, 
respectively. The expressions of J(fcj) in equations (ITS]) . (IT^ and (1201) can be obtained 
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from equation fllll) with the use of the respective Fermi momentum in place of kt. 
The indices /” and ul" are used to denote the interaction between a pair of 
neutrons having the same and opposite spin orientations, respectively. The splitting 
of the strength parameters £g, and of PNM are subject to the condition that 
^o=(^o^ + ^o’'V2, £iy={£’7’‘ + and s^ex={^^Jx + 4^0/2- The exchange strength 


parameter in CPNM in equation fl20l) can be expressed in terms of the hnite range 
strength in the TO state of the N-N interaction as, 


T/2 

^ex PO 0 


f{r)d^r. 


( 22 ) 


The other parameters of CPNM for SEI in equation fl20|) are (due to the zero 

range of the density-dependent term of SEI) and Sq — ~^ex- The SNM is completely 


determined by the parameters 6 , 7 , a and the combinations 


_i_ ^ul 
tQ Lq 


— ^o'l 


p/ _|_ ^ul 

^7 ^7 




7 ’ 


p/ I ^ul 

^ex ' ^ex 


= 


These strength parameters in SNM can also be written as, 

POfo-fiSE , o-fiTE MrTO \rSO\ / 


lb 


0 


V?") / S(r)d\ 


(23) 


(24) 


£0 = \hp, + ^(317'= + Sr™ + 91/™ + V/™) 


f(r)£r, 


(25) 


(26) 


In the foregoing equations, J f{r)d^r — can be replaced where ever it ocures for 

the Gaussian form of f{r). 

We shall now briefly outline the procedure of determination of the parameters, 
as adopted in previous studies of NM and hnite nuclei ET]. The range a and the 
exchange strength in SNM are determined by means of a simultaneous minimization 
procedure using the experimentally extracted constraint [l9l [50], [55] that the attractive 
optical potential changes sign for a kinetic energy 300 MeV of the incident nucleon. The 
NM values of the saturation density po and energy per particle e(po) at saturation are 
the only quantities needed to completely determine a and Sex (see [To] for details). The 
parameter b is hxed for avoiding the supraluminous behaviour in SNM [6T]. It reads 


bpo= 


Mc^ 


Tq/^-^Ro) 


(7+1) 


-1 


with T, 


fPk'i 


fo~- 


fo 


2M 


where kf^ is the Fermi momentum in 


SNM at normal density and M is the nucleonic mass. Its calculation requires again the 
knowledge of the NM values poi e(po) and the parameter 7 . The stiffness parameter 7 
determines the density dependence of the EOS in SNM. The two remaining parameters 
in SNM, namely and are obtained from the satnration conditions, that is, from the 
values of e(po) and po- The stiffness parameter 7 is kept open and its admissible values 
are constrained by the condition that the pressure-density curve lies within the region 
extracted from the analysis of flow data in heavy-ion collisions (HIC) at intermediate 
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and high energies [52] • Thns, a complete stndy of SNM can be performed for a given 7 
if one assnmes standard NM valnes of po and e(po)- 

To extend the stndy for ANM, one needs to know how Sex, and Eq split into like 
and nnlike isospin channels. The splitting of Eex into and is decided nsing the 
physical constraint resnlting from the stndies of the thermal evolntion of NM properties 
[59] . This stndy predicts a critical valne of the splitting of e^x for which the thermal 
evolntion of NM properties as well as the entropy per particle in PNM does not exceed 
that of SNM. The resnlting critical valne of the splitting is = “^eex- The n-p effective 
mass splitting predicted with this choice of nicely coincides with the resnlts of DBHF 
calcnlations | 1 ] as has been shown in the previons work EH. The splitting of the 
remaining two parameters, namely and £95 is obtained by assnming, on the one 
hand, a standard valne of Es{po) at satnration and, on the other hand, the valne of 
its derivative E'g(po) = for which the asymmetric contribntion of the nncleonic 

part of the energy density in charge nentral beta-stable n-|-p-|-e-f p matter (referred to 
as nentron star matter) becomes maximnm. This choice predicts a density dependence 
of the symmetry energy which is neither very stiff nor soft and does not allow the direct 
URCA process to occnr in nentron stars. The popnlation synthesis models [53] based 
on cooling calcnlations [63| predict that there shall be no direct URCA process at least 
in typical nentron stars. Constraining the splitting of the three strength parameters Eex, 
E.y and £0 allows one to determine the nine parameters that describe the ANM. The SEI 
with the parameters obtained in this way has the ability of reprodncing the microscopic 
trends of the density dependence of the EOS and the momentnm dependence of the 
mean helds in ANM [571 EH] • 

There are still two parameters open, which were taken to be to and xq in the 
previons work [37] and were determined from hnite nnclei calcnlations. The energy 
density of a hnite nnclens was constrncted from the nnclear, Conlomb and spin-orbit 
interactions. An improved semi-classical ^^-approximation [65l [ 66 ] was nsed to localize 
the exchange contribntions of the nnclear part. The energy density was thus expressed 
in terms of the local variables, namely, nucleon densities, kinetic energy densities and 
spin-densities. Utilising the variational principle results into Skyrme-like Hartree-Fock 
equations which were solved to get the neutron and proton orbitals. Using this energy 
density functional, directly derived from the SEI, one could determine the two pending 
parameters to and Xo along with Wq, the strength of the spin-orbit interaction, from the 
experimental binding energy of the double-closed-shell nuclei ^°Ca and from 

the lp 3 / 2 -lpi /2 level splitting in as discussed in Ref. [37] . In the calculation of the 
ground-state properties of open-shell nuclei, the pairing correlations were considered 
in the BCS approach with a density dependent zero-range pairing interaction |S7|. 
The binding energies and charge radii of magic nuclei reproduced the corresponding 
experimental values within an accuracy of 0.1%. It was also found that the energy 
density functional associated to the SEI including pairing was well suited for describing 
binding energies and charge radii of open shell nuclei. The experimental binding energies 
of 161 even-even spherical nuclei and the measured charge radii of 86 even-even spherical 
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nuclei from to 224-[j reproduced by our model within windows of ±2 MeV and 
±0.02 fm with overall root mean square {rms) deviations of 1.54 MeV and 0.015 fm, 
respectively. These deviations are in consonance with those obtained in other common 
mean held interactions. For example, as given in Table 3 of Ref. [60], the corresponding 
rms deviations in binding energies and charge radii for the same set of nuclei considered 
here are of 1.71 MeV and 0.024 fm in the Skyrme force SLy4, of 3.58 MeV and 0.020 
fm in the RMF set NL3, and of 2.41 MeV and 0.020 fm in the Gogny force DIS. 

To constrain the two parameters xq and Wq from hnite nuclei using the 
aforementioned protocol, namely, BE of lp 3 / 2 ^ 1 pi /2 level splitting in 

shows some arbitrariness in the sense that small variations in Wq that imply an 
appropiate change in xq, do not modify the quality of the overall rms deviations of 
binding energies and radii. For example, xo=0.6 and Wo=llb MeV give rms deviations 
of binding energies and radii of 1.54 MeV and 0.015 fm, whereas the combination 
Xo=0.2 and Wo=116 MeV gives rms deviations of 1.47 MeV and 0.015 fm. It is to 
be noted here that to determined from the BE of “^^Ca is almost insensitive to the 
choice of Xq (because N = Z) and Wq (because ^°Ca is spin saturated). However, 
the uncertainty in determining xq from hnite nuclei has an important impact on the 
individual contributions to the SE, TE, TO and SO states of iV — iV interaction in 
NM. This uncertainty largely manifests in spin properties of NM, as shall be discussed 
in the forthcoming section. 

3. Results and Discussions 

It is evident from the discussions of the foregoing section, on the procedure of hxing the 
parameters of SEI, that for a given value of 7 a complete study of SNM, PNM, ANM 
and hnite nuclei can be performed only by assumming standard values of po, e(po) and 
Es{po). The empirical values of these three NM properties used by diherent models vary 
over certain ranges as mentioned in section 1 . Out of these three properties, e(po) has 
the minimum uncertainty and for all the models its value lies in the range —16 ± 0.2 
MeV. As regards the other two NM properties, the majority of the RMF sets have a 
value of the saturation density po to the lower side of the range 0.17 ± 0.03 fm“^ and a 
value of the symmetry energy Es{pq) in the higher side of the range 33±5 MeV. In the 
non-relativistic mean held theories the value of po centers around 0.16 fm“^ and Es{po) 
to the lower side of the above mentioned range. In both relativistic and non-relativistic 
microscopic calculations po is predicted in the higher side of its range. 

In the present work we shall see that po is strongly correlated with 7 that determines 
the stihness of the EOS in NM. In this context we have hrst examined the pressure- 
density {P ~ p) relation in SNM for diherent values of 7 . The pressure in SNM is 
calculated as P{p) = from equation (] 8 |), where we have used the values e(po)=—16 

fPiP 

MeV and Fermi kinetic energy Tf^= MeV (corresponding to po=0.161 fm“^), 

with kfQ = { being the Fermi momentum at saturation density. The P ^ p curves 

for 7=1/6, 1/3, 1/2 and 2/3 corresponding to incompressibility of NM, iF(po)=207, 226, 






J. Phys. G: Nucl. Part. Phys. 


T. R. Routray et al. 10 


Table 2. Values of nuclear matter properties at saturation (namely, incompressibility, 
Fermi kinetic energy, density and symmetry energy), rms deviations from experiment 
in energies and charge radii, and the parameters to, xq and Wq for the four EOSs of 
SEI corresponding to 7 = 1 / 6 , 1/3, 1/2 and 2/3. 


7 

A'(Po) 

MeV 

Tfo 

MeV 

Po 

fm“^ 

Es{po) 

MeV 

SE 

^-‘-'rms 

MeV 

^^rms 

fm 

to 

MeVfm^ 

Xo 

hho 

MeV 

1/6 

207 

37.2 

0.1623 

36 

1.4916 

0.0187 

-575 

-0.7 

120 

1/3 

226 

36.8 

0.1597 

35.5 

1.3932 

0.0167 

201 

1.1 

118 

1/2 

245 

36.4 

0.1571 

35 

1.5402 

0.0152 

437 

0.6 

115 

2/3 

263 

36.1 

0.1552 

35 

1.9336 

0.0154 

540 

1.38 

112 


245 and 263 MeV, respectively, are verified to pass within the experimentally extracted 
region of Ref. j62]. We shall now obtain all the nine parameters of ANM for each EOS 
of these fonr 7 by assumming Es{pq)=33 MeV together with the values e(po)=-16 MeV 
and TfQ=37 MeV used above for SNM. Then the study of finite nuclei is performed for 
each of these four EOSs by adopting the procedure for the determination of R, xq and 
Wq outlined in the last section. The results for the deviations with respect to experiment 
in the binding energies, 6E, of 161 even-even spherical nuclei and in the charge radii, 
Svch, of 86 even-even spherical nuclei are calculated for the four EOSs corresponding to 
7 = 1 / 6 , 1/3, 1/2 and 2/3. The results of 6rch of the four EOSs presented in Figures 
1. The rms deviations from experiment in the charge radii 6rrms in hgure 1 for the 
four EOSs corresponding to 7 = 1 / 6 , 1/3, 1/2 and 2/3 are 0.036 fm“^, 0.017 fm“^, 0.038 
fm“^ and 0.057 fm“^, respectively. From this figure, a strong correlation between 7 and 
Po on the radius of finite nucleus can be seen. It can be concluded that for a given 
stiffness of NM there is a critical value of po for which the deviations in the charge 
radii of all these nuclei center around zero giving a minimum rms value, Sr rms, of the 
deviations. A relatively softer (stiffer) EOS with po corresponding to Tf^=37 MeV 
overestimates (underestimate) the radii in finite nuclei. It has been verified that this 
conclusion does not change on the choice of either the e(po) or the Eg^po) values. With 
the four considered EOSs we have made several calculations by varying e(po) and Es{po) 
and have searched for the optimal rms results for the BE and radii of the considered 
set of spherical nuclei, following the procedure outlined in section 2. The results of the 
minimum rms deviations in BE, SErms, and in charge radii, Sr rms, some nuclear matter 
properties as well as the SEI parameters to, Xq and Wq for the four considered EOSs 
are given in Table 2. The remaining parameters of the SEI corresponding to these four 
EOSs are given in Table 3. 

The results of po in table 2 reveal that for a softer EOS, the central density is 
predicted to have a relatively higher value as compared to that of a stiffer EOS. A similar 
behaviour is also observed for the symmetry energy. As the incompressibility of NM 
changes from 207 MeV to 263 MeV, the saturation density decreases from 0.1623 fm“^ to 
0.1552 fm“^ (corresponding to a decrease of the Fermi kinetic energy from 37.2 MeV to 
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Figure 1. Deviation in charge radii, 5rch, for 86 even-even spherical nuclei with 
nucleon number between A = 16 and A = 224 for the four EOSs corresponding to the 
same values of Tf„=37 MeV and Egipo) =33 MeV but different 7=1/6, 1/3, 1/2 and 
2/3 in (a), (b), (c) and (d), respectively. 


Table 3. Values of the parameters of asymmetric nuclear matter for the four EOSs. 


7 

b 

fm^ 

a 

fm 

^ex 

MeV 

MeV 

^0 

MeV 

MeV 

MeV 

6^ 

MeV 

1/6 

0.2720 

0.7568 

-96.8427 

-64.5618 

-215.8954 

-131.0415 

213.3364 

142.7149 

1/3 

0.4184 

0.7582 

-95.6480 

-63.7653 

-112.7493 

-67.0819 

110.7436 

78.7768 

1/2 

0.5914 

0.7597 

-94.4614 

-62.9743 

-78.7832 

-45.8788 

77.5068 

57.7687 

2/3 

0.7852 

0.7609 

-93.5766 

-62.3844 

-61.9929 

-33.9536 

61.6896 

47.0768 


36.1 MeV), and the symmetry energy changes from 36 MeV to 35 MeV. This correlation 
between K{pq) and po is also found in earlier Skyrme II-VI sets [25] and conforms to 
the fact that as the matter becomes stiffer, the internucleon separation increases. This 
conclusion is also substantiated if one examines the values of the incompressibility and 
the saturation density of various parameter sets of the RMF model that are successfully 
applied to hnite nuclei calculations (e.g., the popular NL3 set [IH] has a saturation 
density 0.148 fm“^ and an incompressibility 271 MeV). 
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Figure 2. (Color online) {a){upper pannet) The neutron-proton effective mass splitting 
in neutron-rich ANM shown as a function of the isospin asymmetry /3^ at normal 
density for the four EOSs of table 2. ih){lower pannet) The density dependence of the 
symmetry energy is given for these four EOSs. 


The variations in the valnes of po and Es{po) with the stiffness of NM as ascertained 
from the stndy of hnite nnclei have rather small inflnence on the ANM resnlts. This is 
shown in Fignre 2, where the n-p effective mass splitting {m*/m)n-p in ANM at normal 
density as a fnnction of the isospin asymmetry, (3r ={pn — Pp) / (Pn + Pp), and the density 
dependence of symmetry energy are given in the npper and lower panels, respectively, 
for the fonr sets of EOSs. The identical resnlts for the n-p effective mass splitting for 
all fonr EOSs can be nnderstood from the almost same valnes of the exchange strength 
parameters and range a, given in table 3, for the EOSs. These parameters 

determine the momentnm dependence of the mean helds in SNM and PNM. In the 
process of determination of the exchange strengths Sex and in SNM and PNM, 
two particnlar combinations and (Vq^ — 31/),™) of the strenghts of the 

hnite range part of the interaction in the fonr basic states of the N-N interaction are 
getting hxed (see eqnations ([6]) and ([7])). These two combinations in terms of the known 
interaction parameters are given as, 

fo-{rTE ■{rSO\ _ 


(27) 
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SSex 

3(po7r^/^Q!^) ’ 


(28) 


where we have used that = ^Sex and = ^Sex- 

The four EOSs considered here having different stiffness have nearly identical values 
of Eex (and ^gj.) and, therefore, shall have a similar momentum dependence of the mean 
helds in SNM (and PNM). The small differences are due to the variation in their values 
of po which he within a close range. Therefore the combinations (3 Vq^® — Vq^) and 
{Vq^ — ?>Vq^) have similar values for all the four EOSs. 

In the course of the determination of the strength parameters eo in SNM and 
£q in PNM, the interaction parameters involved are mutually adjusted subject to the 
constraint that {?)Vq^ — Vq^) and (Vq^^ —3V(,™) retain the values specihed by equations 
(l27|l and (l28ll . This results into unconstrained variations of the strengths in the four 
states of the N-N interaction for the four different EOSs considered. This is shown in 
Figure 3 where the individual contributions to the interaction part of energy per particle 
in SNM, < V > /A, m the four basic states of iV — iV interaction are shown as functions 
of the Fermi momentum. It can be seen that the contributions do not follow any dehnite 
pattern with respect to their NM properties, in particular, with increase in stiffness of 
SNM. The TO state contribution of the EOS with 7 = 1/3 is predicted to be attractive, 
whereas it is repulsive for the three other EOSs, being more attractive for 7 = 1/6 than 
in 1/2 case. Similar arbitrary behaviour is observed in SO and TE channels, where, the 
curves of the EOSs of different stiffness do not follow any dehnite trend. 

As mentioned in the previous section, for a given EOS, a small readjustment in 
the values of the parameters Xq and Wq (subject to the constraints of the BE of ^ospp 
and the lp 3 / 2 ^ 1 pi /2 level splitting in ^®0) allows to describe BEs and charge radii of 
spherical nuclei with similar rms deviations. For example, a:o=0.2 and Wo=llQ MeV 
could be an alternative set for the one given in table 2 (a;o=0.6 and Wq=115 MeV) for 
7 = 1 / 2 . Though these two EOSs for the same 7 give identical results in the isospin 
channel of ANM (same n-p effective mass splitting and same density dependence of the 
symmetry energy), the behaviour of their contributions in the four basic channels of 
the N-N interaction are found to differ apreciably. This is manifested particularly in 
the spin channel and the spin symmetry energy Eo-(p) (calculated using the expression 
in equation (17) of our earlier work [37]). This is illustrated in Figure 4, where the 
contribution of the TO state to the interaction energy in SNM and the spin symmetry 
energy are plotted as functions of the density in the upper and lower panels, respectively, 
for the two sets of values for xq and Wq corresponding to the EOS of 7 = 1 / 2 . The 
same situation happens for each of the four EOSs considered in the work. It may be 
noted that the predictions in the spin channel are crucial in the studies of magnetic 
properties of dense NM. Spin polarisation properties in various types of NM have been 
studied extensively in theoretical approaches using both microscopic and effective models 
(MlESlEQliniESlEslIHlZSllTniEZllZHllZnilHnilHIllHaEsllHlI often with contradictory 
conclusions. The spin polarizability of NM can have strong effects on the neutrino mean 
free path and can impact on the formation mechanism and cooling scenario of neutron 
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Figure 3. (Color online) Contributions to the interaction part of the energy per 
particle, < V > jA^ in SNM coming from the four basic states of the N-N interaction, 
SE^ TE, TO and SO, as a function of the Fermi momentum for the four EOSs 
corresponding to 7 = 1 / 6 , 1/3, 1/2 and 2/3. 


stars. 

The divergent predictions in the spin channel of the different EOSs for a given 7 
value arises due to the freedom allowed in the determination of the two parameters xo 
and Wq from hnite nuclei without compromising much the predictions of BEs and radii. 
This results into several possible values of Xq, and hence different sets of values for 

and subject to the constraint that the combinations (3 Eq^® — and 
— ?>Vq^) remain invariant. A possible way to remove the existing arbitrariness is 
to determine the parameter xq from NM keeping only Wq open for hnite nuclei (apart 
from to that is hxed from BE of ^°Ca) provided the overall predictions in hnite nuclei 
do not worsen. 


3.1. Determination of xo 

The parameter xq can be expressed as, 


Poto 


(^0 “ 4a; + 0^4); 


xo = l- 


(29) 
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Figure 4. (a){upper panel) Contributions to the interaction part of the energy per 
particle < V JA in SNM coming from the TO state as a function of the Fermi 
momentum for the two EOSs having the same 7=1/2 but differing in their xq and Wq 
values, {h){lower panel) The density dependence of the spin symmetry energy for the 
same two EOSs of 7=1/2. 


from equations o. o and ([22]). With R fixed from finite nucleus and Eq and 
known from ANM, we can determine xo with the knowledge of The parameter 
decides the momentum dependence of the mean held in CPNM and can be ascertained 
from the splitting of the parameter of PNM between spinwise like and unlike 
pairs in PPNM. This splitting is subject to the constraint 25^^. Thus the 

splitting into the spinwise-like channel, can take any value between 0 and 
and correspondingly is decided. We have examined the effective mass splitting, 
{m*/m)nu-nd = {tTP/'rn)nu — {^*/fri)nd, between spin-up (nu) and spin-down (nd) 
neutrons in PPNM at normal density po for various possible values of and have 
compared the results with the DBHF prediction with the Bonn B potential [75]. The 
effective masses of nu and nd neutrons in PPNM can be calculated from the expression, 

771 dUnu,nd{}^i Pi 


Tfl 

m 


nu.nd 


1 + 


n^k 


dk 


(30) 
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where Unu(nd) is the mean field of nu {nd) neutrons in PPNM and Pu is the spin 
asymmetry defined as (3„={pnu — Pnd)/Pi with p^u and pnd being the densities of nu 
and nd neutrons and p = pnu + Pndi the total density of PPNM. The mean fields of nu 
and nd neutrons in PPNM calculated with the SEI give the results, 


and 


/ 7 O \ _ ^0 I ^0 I ^7 f ^ \ ^ 

'^nuy^ 1 P 1 Pa) pnu T Pnd T ~~piI\Pnd \ 1 H | 7 

Po Po pI ^ P I+ hp 

p/,/ ^l,ul 

H- —Pnul{ki knu) H- —Pndl{k, knd) 

Po Po 



7 


(31) 


Pi / 3 (t) 


A,I Lul A.ul 

_ ^0 ^0 ^7 

— - Pnd H- P nu + nu (IT 

Po Po p^^ V 


pnd 7 

p 1 + hp) 


P,l 


-Lul 


H- —Pndl{ki knd) H —Pnul{ki knu)i 



7 


(32) 


Po Po 

where I{k,ki) {ki = k^u, knd) is given by equation ( 1 ^ and it has been taken into 
account that = Q as discussed in Section 2. Using the SEI mean fields in equation 
(130]) it is found that the nu effective mass, {m*/m)nui becomes larger than the nd 
effective mass, {m*/m)ndi for within 0 and which is the trend observed in the 
microscopic DBHF calculation. The difference between the nu and nd effective masses 
is maximum for = 0 and this difference decreases and becomes zero as increases 
in magnitude and coincides with Beyond this value and up to = 2e[^, the 
trend reverses and the nd effective mass becomes larger than the nu effective mass. 
Upon comparison with the DBHF result, it is found that for the particular value of the 
splitting £T=£g 3,/3 the results of {m*/m)nu-nd for all the four EOSs, corresponding to 
the four 7 values, are in close agreement with the microscopic DBHF prediction ua over 
a wide range of the spin asymmetry This is shown in Figure 5, where the calculated 
results of {m*/m)nu-nd at saturation density for the value e’‘^lj.=e[^/3 along with the 
DBHF prediction are plotted as a function of the spin asymmetry /3„. The momentum 
dependence of the mean field in spin polarized PNM is fixed once is known. This, in 
turn, provides us the finite range strength in the triplet-odd state which is obtained 
from equation (|2^ using the given value = e[^/3. The V/f^ strength resulting from 
this splitting, i.e., Uq™ = —£ex/(P ott^'^^o^) is found to be repulsive is always negative 
in table 3), and ranges within V/f^ = 54.9 — 54.6 MeV as 7 varies from 1/6 to 2/3. The 
repulsive character of V/f^ is an essential requirement for the effective mass in CPNM 
at normal density to be smaller than 1, as well as for the stability of CPNM at any 
higher density. 

The parameter xq, which was determined from the binding energy of in our 

earlier fitting procedure [37], now connects both NM and the finite nucleus. In the 
present work, once we know from NM and to from the finite nucleus, the parameter 
Xq is determined from equation fl29l) . We fix to from the BE of ^°Ca and the spin-orbit 
strength Wq from the BE of ^°®Pb. The resulting parameters are reported in table 4. 
The to values corresponding to different stiffness 7 of the EOS remain the same as in 
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Figure 5. (Color online) Effective mass splitting, {m* /m)nu-nd, between spin-up 
and spin-down neutrons in polarised PNM as a function of the spin asymmetry 
Pa={pnu — Pnd)Rpnu + Pnd) for the four EOSs with eC=eg3,/3. The result of the 
microscopic DBHF calculation m is also given for comparison. 


table 2 , whereas Wq changes slightly owing to the change of the xq parameter, bnt 
being free from any arbitrariness. The contribntions to the interaction energy in SNM 
in the fonr channels resnlting from this parametrisation are shown in Fignre 6 . Upon 
comparison with the earlier resnlts shown in Fignre 3, it can be seen that a systematic 
variation with respect to the stiffness of the EOS in SNM is obtained in the present way 
of determining the parameters. The interactions in the TO and TE channels now have a 
similar behavionr for all the fonr EOSs, whereas in the SO and SE channels systematic 
variations with respect to the stiffness parameter 7 are observed. In both the SO and 
SE states, the EOS having lower incompressibility predicts relatively more attraction. 
In the same fignre the resnlts obtained with different Gogny parameter sets are given 
for comparison. With the new method of hxing R and Xq for SEI, the arbitrariness in 
the prediction of the spin symmetry energy E^{p) (see lower panel of hgnre 4) for each 
EOS also gets removed. The finite range strengths in the fonr basic channels Vq^, 
and can be expressed, with the choice £g’,^=£g 2 ,/ 3 , in terms of R as follows: 


{pQTT^Oa^) ’ 


( 33 ) 
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Figure 6. (Color online) Contributions from the four basic states of N-N interaction 
SE^ TE, TO and SO to the interaction energy in SNM as a function of the Fermi 
momentum for the four EOSs corresponding to 7 = 1 / 6 , 1/3, 1/2 and 2/3 having the 
parameter Xq determined from the study of spin polarised PNM and given in table 4. 
The results of the Gogny Dl, DIS and DIM forces are also shown. 




yTE 

*^0 


= 


2e. 


ex + 8 ^ 0 ) ~ StoPo 

(lO^ea; + S^o) ~ S^oPo 
(poTT^/^Q!^) 


(34) 

(35) 

(36) 


Once all the parameters of SEI have been hxed, the mean helds Unu{k, p, /3a) and 
Undik, p, /3a) of uu and nd neutrons in PPNM can be calculated as a function of the 
momentum k. The mean helds for the SEI set of 7 = 1/2 are shown in Figure 7(a) at 
density po and spin asymmetries /3a=^, 0.1, 0.2 and 0.4. The curves for nu and nd 
neutrons lie above and below the unpolarized curve almost symmetrically, with larger 
separation between them for larger /3a. The dependence of the Unu and Und mean helds 
on the spin asymmetry, calculated at p=po and k=kn for the four EOSs corresponding 
to 7 = 1 / 6 , 1/2, 1/3 and 2/3, is shown in Figure 7(b) along with the DBHF prediction 
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Table 4. Values of the parameters xq and Wq, where xq is determined from the 
considerations of both NM and finite nucleus, given in equation (1^^ . for the four EOSs 
along with their rms deviations in BE of 161 even-even nuclei and charge radii of 86 
even-even nuclei. 


7 

^0 

MeVfm^ 

Xo 

VPo 

MeV 

^^rms 

MeV 

A'T 

rms 

fm 

1/6 

-575 

-1.1161 

118 

1.6993 

0.0189 

1/3 

201 

3.1928 

115 

1.6754 

0.0170 

1/2 

437 

1.4192 

112 

1.8518 

0.0155 

2/3 

540 

1.0659 

113.5 

1.8297 

0.0152 


m- In order to have a direct comparison of the results, the curves have been shifted 
to the origin by subtracing the respective Unu(nd){kn, Po, (3a = 0) values. The behaviour 
of the mean helds of all the four EOSs of SEI is alike but their splittings are smaller in 
comparison to the DBHF prediction. The SEI result, with the value and 

the DBHF result have a closely similar momentum dependence in their mean helds in 
PPNM, as evident from the effective mass property in hgure 5. Thus, the difference of 
the results in hgure 7(b) can be attributed to the density-dependent part of the mean 
helds, that is largely accounted for by the spin symmetry energy. In the case of SEI, the 
splitting of the strength parameter of PNM, corresponding to the density-dependent 
part of the interaction, into the like-spin channel vanishes, i.e., = 0, due to the 

zero range of the density-dependent term. As a consequence, the energy per particle in 
CPNM cannot have a stih enough behaviour to ensure that the neutrons in polarized 
state shall have higher value at all densities as compared to unpolarized state, the trend 
obtained in the microscopic BHF and DBHF calculations [73l[75]. This indicates that 
the density-dependent part of the interaction in the case of SEI needs to be improved in 
order to reproduce with better quality the microscopic trend of the density-dependent 
contribution in spin polarized matter, which shall not be considered in the present work. 


In order to have further insight into the momentum dependence of the mean helds 
in PPNM and the effective mass splitting, we note that the sign of the effective mass 
splitting {m*/m)nu-nd is decided by the dimensionless quantity 


m p, (3a) 

fpk dk 


(37) 


where p, f3a) is the difference between the nu and nd mean helds given by 

equations m and ([32D- If equation fimi is negative (positive) the nu (nd) ehective 
mass is larger than the nd (nu) ehective mass, as can be easily deduced from equation 
In the limit (3a 0, equation (15711 reads as 


m '^’^{k,p,(3a 


rpk 


dk 


k — kn 


= -(3c 


^ex ) 


pA^ m 


1-1 + 


4^ 

A2 


4fc 

e~^ 


( 38 ) 
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Figure 7. (a) Mean fields of nu and nd neutrons in PPNM as a function of momentum 
k at density pQ for the EOS corresponding to 7=1/2 having values of spin asymmetry 
/3 ct= 0, 0.1, 0.2 and 0.4. (b) Mean fields of nu and nd neutrons in PPNM at density 
Pq and momentum k=kn as a function of spin asymmetry (3^ for the four EOSs 
corresponding to 7=1/6, 1/3, 1/2 and 2/3 alongwith the DBHF results [75] . 


This equation can predict the nature of the effective mass splitting in spin polarized 
PNM by inspecting the value of Due to the fact that the square bracketed 

factor in equation (138|) is positive and the exchange strength parameters are attractive, 
the nu {nd) effective mass will lie above the nd {nu) effective mass if {s^ex~^^ex) D positive 
(negative). In terms of the interaction parameters, we have (^e^— 

and hence the sign of (j38jl can be predicted in terms of the value oiW — H. Therefore, 
to have a nu effective mass larger than the nd effective mass requires that W — H he 
negative. In the case of the Gogny interaction, the short-range term being dominant, 
the behaviour of the nu and nd effective masses can be predicted from the value of 
Wi — Hi (i.e., the quantity hh — if in the short-range term of the Gogny forces). In the 
Gogny Dl, DIS, DIN and DIM parameter sets, Wi — Hi assumes positive values and 
it has been verihed that in these sets the nd neutron effective mass in PPNM lies above 
the nu effective mass. 
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3.2. Spin polarized SNM 


We shall study the two extreme cases of spin polarisation, CSNM and CASNM, in 
SNM. The energy per particle in CSNM and CASNM using the SEI is calculated from 
equations ([9]) and ffTOj) and is shown in the Figures 8. (a) and (b), respectively, for the 
four EOSs corresponding to 7=1/6, 1/3, 1/2 and 2/3 together with the results of the 
SNM. The results in the hgures are shown upto 10 times the normal NM density po as 
the central densities of the maximum mass neutron stars obtained in this range of NM 
incompressibility are found, in an earlier work 1571, to entend upto 8 to 10 times the 
normal density. Like the isospin symmetry energy in ANM, the spin symmetric energy 
in these two types of spin polarized SNM can be expressed as the differences of energy 
per particle in SNM from that of in CSNM (CASNM), 


E'M = - '(p). (39) 

with i=S(AS). The spin symmetry energy thus calculated for these two cases of 
polarisations in SNM are shown in Figures 9. (a) and (b) for the four EOSs of SEI 
alongwith the results of microscopic and effective models. The spin symmetry energy for 
SEI calculated from equation (17) of Ref. [37] . derived under the Taylor series expansion 
of the energy per particle in spin asymmetric ANM, is also shown in hgure 9. (a) (curves 
with crosses) and is seen to compare well with the results of CSNM. 

It can seen from hgure 8. (a) and 9.(a) that the SEI force sets do not predict EM spin 
ordering in SNM in agreement with the predictions of all microscopic models |H2l [77] 
and some of the effective models. However, AFM spin ordering in SNM is predicted 
by the SEI sets as in case of Gogny DIS force (cf. Ref. [79]), which is contrary to the 
microscopic predictions. The SEI sets allow the antiferromagnetic transition at a density 
about four times the saturation density in SNM. The results for the three considered 
Skyrme sets, namely SLy4, SLy9 ans SkI3 [26] ET] show divergent behaviour. For the 
SLy4 and SLy9 forces, a EM transition is not realized at any density, whereas AFM 
ordering of the neutron and proton spins in SNM occurs at relatively high densities. 
In the case of SkI3 both EM and AFM spin ordering are possible, with the former 
being realized at a relatively smaller density than the latter. The Skyrme results are 
calculated from equation (13^ . where the expressions of the energy per particle in CSNM 
and CASNM are. 
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Figure 8. (Color online) (a) Energy per particle in CSNM in equation @ is shown 
as functions of density p for the four EOSs of SEI corresponding to 7 = 1 / 6 , 1/3, 1/2 
and 2/3, with the sets of parameters determined for the values of Sq in table 4. (b) 
Same as (a) but for CASNM in equation (nni). 


(41) 


respectively. 

For the sake of further insight, we express the energy per particle of CSNM and 
CASNM given in equations (jO]) and ffTOj) in terms of the hnite range strength parameter 
in the triplet-odd state Vq^ together with the other known parameters, as given by. 
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Figure 9. (Color online) (a) Spin symmetry energy, E^{p), calculated from 
equation (15^ for FM spin polarised SNM as a function of density for the four EOSs 
corresponding to 7=1/6, 1/3, 1/2 and 2/3 with the sets of parameters determined 
for values of xq in table 4. The corresponding results calculated from the equation 
(17) of ref. [37] are shown by curves with crosses. The results of the microscopic LOCV 
calculation with AVis interaction [53] and of the Gogny Dl, DIS and DIM [5011371135] : 
Skyrme SLy4, SLy9 and Skl3 force sets [laiiT] are also shown, (b) Same as (a) but 
for the AFM spin polarised SNM. 
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respectively, where can be obtained from equation ([TT|) for ki=kK^\ Now, we can 
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write them as 

(44) 

where 44 s(as)(p) and Bs{as){,p) represent the hrst and second sqnare bracketed terms in 
eqnation (I42|) (eqnation (l43D i. repectively. The qnantities and Bs(as) are shown 

in Fignres 10.(a) and (b), respectively, as a fnnction of density for both CSNM and 
CASNM. The results in these two hgures show that the energy per particle in CSNM, 
epoiip), shall be an increasing function of density for repulsive The procedure 

adopted in the determination of xq, discussed in the foregoing subsection, also predicts 
a repulsive value of and hence an increasing trend of the spin symmetry energy is 
obtained for the EOSs of SEI, exhibiting stability against EM spin ordering in SNM. The 
results of As{p) for the four EOSs obtained under the present procedure of determination 
of the parameters are almost identical. Hence, the CSNM energy per particle epoiip) 
shows little dependence on the stiffness of SNM, as can be seen from the curves for the 
four EOSs of CSNM in hgures 8 .(a) and 10.(a). Due to the fact that e(p) is stiffer for 
higher 7 values, the spin symmetry energy Efj{p) shall exhibit a softer behaviour and 
this can be seen from hgure 9(a). On the other hand, the energy per particle in CASNM, 
epj(p), has a dependence on the stiffness parameter 7 of the EOS, but this dependence 
is not stronger than its counterpart in SNM as one can realize from hgure 8 .(b). From 
hgure 10 .(a) it can be observed that Aas{p) has a stiher behaviour for higher 7 values 
but it remains below its counterpart curve As{p) of CSNM at all densities. Moreover, 
the contribution of the Bas{p) term is negative for a repulsive as can be seen from 
equation (1T3|1 and hgure 10. (b). This makes the energy per particle e^i(p) in CASNM 
softer than its counterpart e^^i^p) in SNM, and AFM spin ordering is predicted at a 
critical density close to 0.65 fm“^ for all the four EOSs of SEI, as found in hgure 9.(b). 
It may be pointed out that the AFM spin polarized SNM can be stable against the 
AFM transition for an attractive Vq^. But in that case, the FM spin polarized SNM 
shall be realized. It may be mentioned that FM and AFM spin ordering in diherent 
types of NM is still under debate and predictions of diherent model calculations are 
often contradictory (cf. Ref [83] ). 

3.3. Finite nuclei properties with the new parameter sets 

At last we examine the ability of SEI in predicting hnite nuclei properties with the new 
sets of parameters where Xq is hxed from the consideration of spin polarized neutron 
matter as explained in section 3.1. In this new protocol, the spin-orbit strength Wq 
is adjusted to reproduce the BE of ^°®Pb. The prediction of the deviations in BEs 
and charge radii for the even-even spherical nuclei are shown in Figures 11 and 12, 
respectively, for the four EOSs. The rms deviations in BE, SE^ms, and charge radii, 
hrrms, along with the values of R, xq and Wq are listed in table 4. From Figures 11 and 
12 and the corresponding rms deviations given in Table 4, it can be seen that the results 
of hnite nuclei are reproduced within reasonable accuracy for all the four EOSs. The 
present values of the SErms and 6rrms deviations are comparable to those obtained for 
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Figure 10. (Color online) {(a) upper panel) The function ^s(as)(p) defined in equation 
(1421) (equation (l43l) l for CSNM (CASNM) is shown as function of density p for the four 
sets of EOSs corresponding to 7=1/6, 1/3, 1/2 and 2/3 with the sets of parameters 
determined for values of xq in table 4. ((b)lower panel) The function -Bs(as)(p) defined 
in equation (1^ (equation (l43l) l for CSNM (CASNM) as function of density p for the 
same four EOSs as in (a). 


the same set of nuclei with the SLy4, NL3, DIS, BCPl and BCP2 effective forces given 
in Table 3 of Ref. [60]. In order to examine the predictions of single-particle levels and 
their splittings, the neutron and proton energy levels in ^ospp shown in Figure 13 
for the four EOSs having 7 = 1 / 6 , 1/3, 1/2 and 2/3 alongwith the experimental values 
[HSllSS]- The quality of the single-particle spectra obtained in the present case is similar 
to that obtained in other traditional effective forces shown in our earlier work [37]. From 
the comparison of the spectra in the four EOSs, it can be concluded that nearly similar 
spectra are predicted by all the four EOSs, having a tendency of widening the gaps 
between the single-particle levels with increase in stiffness. 
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Figure 11. Deviation in energy for 161 even-even spherical nuclei with nucleon number 
between A = 16 and A = 224 for the EOSs corresponding to 7=1/6, 1/3, 1/2 and 2/3 
with sets of parameters for xq in table 4 and are shown in (a), (b), (c) and (d), 
respectively. 


4. Summary and Conclusions 

In this work we have investigated some of the various types of nuclear matter as well as 
hnite nuclei with the hnite range simple effective interaction. This interaction describes 
the N-N force in the four basic states SE, TE, TO, SO with the same range but different 
strengths. The SEI contains in total twelve parameters. Most of these parameters are 
determined from the consideration of the basic properties of asymmetric nuclear matter 
and the momentum dependence of the mean held. In previous work [37] , the remaining 
parameters—namely, to, xq and the spin-orbit strength Wo —were htted to a few magic 
nuclei. However, on examining the predictions of the SEI in spin polarized symmetric 
nuclear matter and pure neutron matter, it is found that xo and Wq are correlated in the 
earlier htting procedure that leaves some uncertainty in the values of these parameters. 
This implies that changes in the xq parameter can be compensated by small changes 
in Wo that practically do not affect the rms deviations from the experimental binding 
energies and charge radii of spherical nuclei, but manifests in the predictions in spin 
polarized matter (see e.g. hgure 4). To overcome this difficulty we connected the 
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A 

Figure 12. Deviation in charge radii for 86 even-even spherical nuclei with nucleon 
number between A = 16 and A = 224 for the EOSs corresponding to 7=1/6, 1/3, 1/2 
and 2/3 with sets of parameters for Xq in table 4 and are shown in (a), (b), (c) and 
(d), respectively. 


determination of xq to the microscopic trend of the momentum dependence of the mean 
held in spin polarized PNM. The assessment of Xq was decided from the DBHF nu- 
nd effective mass splitting in polarized PNM [75] and to, determined from the binding 
energy of ^°Ca, and hence acting as a connection between nuclear matter and the hnite 
nucleus. Finally, the spin-orbit strength Wq was htted to the binding energy of ^°®Pb. 
The parameter sets of SEI thus obtained are able to reproduce the hnite nuclei results 
with a similar quality to other traditional forces along with dehnite predictions in isospin 
and spin polarized matter. 

The momentum dependence of the spin-up and spin-down neutrons and the ehective 
mas splitting in spin polarized PNM is discussed. It is shown that the spin-up neutrons 
have larger ehective mass than the spin-down neutrons when the strength of the 
exchange part in the like channel is smaller in magnitude than in the unlike channel, 
and vice-versa. In terms of the interaction parameters, this can be interpreted from 
the fact that if the diherence between the Wigner and Heisenberg strength parameters 
is attractive (repulsive) then the single-particle potential and the ehective mass for 
spin-up neutrons in PPNM shall be above (below) the same quantities for spin-down 
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Figure 13. The single particle levels for neutrons and protons in for the EOSs 

corresponding to 7=1/6, 1/3, 1/2 and 2/3 with sets of parameters for xq in table 4 
The experimental results are also given for comparison. 


neutrons. The SEI is found to be competent to reproduce the microscopic trends of the 
momentum dependence in NM, but requires modihcation in the density-dependent part 
for predicting the spin symmetry energy in closer agreement with the microscopic result. 
With the present density dependence of SEI the ferromagnetic and anti-ferromagnetic 
cases of spin polarized SNM have been examined. It is found that the SEI does not 
predict a ferromagnetic transition of SNM, whereas an anti-ferromagnetic transition 
is allowed at a density around four times the normal saturation density. This result 
is qualitatively similar to the result predicted by the Gogny interaction in the Fermi 
liquid formalism [79]. Similar results are also observed in the case of SLy Skyrme force 
sets. However, the SkI3 Skyrme set allows both ferromagnetic and anti-ferromagnetic 
transitions in SNM, but favouring the former at a relatively small density in comparison 
to the latter. The predictions of SEI are analysed analytically and it is found that 
stability against a ferromagnetic transition requires the interaction in the TO state 
to be repulsive, opposite to the requirement of an attractive TO state interaction for 
stability against an anti-ferromagnetic transition. 

The simultaneous study of nuclear matter of different types and hnite nuclei with 
the SEI practically requires standard values of the three NM properties po, e(po) and 
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Eg{po). From the global description of energies and radii of finite nuclei with different 
sets of EOSs of SEI, the influence of the incompressibility K{pq) stands out in terms of 
the variations in the values of and Es{po) given in table 2. As the incompressibility 
K{pq) varies from 207 MeV to 263 MeV, the Fermi kinetic energy Tf^ decreases from 
37.2 MeV to 36.1 MeV and the symmetry energy Es{po) decreases from 36 MeV to 
35 MeV in order to reproduce the results of BEs and charge radii with minimun rms 
deviations. Thus, the SEI with the parametrization determined in this work can be used 
in the study of asymmetric NM and spin polarized NM as well as in hnite nuclei. We 
have restricted our analysis of hnite nuclei to even-even spherical systems. It is therefore 
important to extend this set of nuclei by including deformed nuclei. Work to adapt our 
numerical codes for deformed calculations is being undertaken and will be the subject 
of future communications. 
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